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Abstract 

We investigate the renormalization properties of minimally doubled fermions, at one loop 
in perturbation theory. Our study is based on the two particular realizations of Borigi-Creutz 
and Karsten-Wilczek. A common feature of both formulations is the breaking of hyper-cubic 
symmetry, which requires that the lattice actions are supplemented by suitable counterterms. 
We show that three counterterms are required in each case and determine their coefficients to 
one loop in perturbation theory. For both actions we compute the vacuum polarization of the 
gluon. It is shown that no power divergences appear and that all contributions which arise 
from the breaking of Lorentz symmetry are cancelled by the counterterms. We also derive 
the conserved vector and axial- vector currents for Karsten-Wilczek fermions. Like in the case 
of the previously studied Borigi-Creutz action, one obtains simple expressions, involving only 
nearest-neighbour sites. We suggest methods how to fix the coefficients of the counterterms 
non-perturbatively and discuss the implications of our findings for practical simulations. 
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1 Introduction 



Minimally doubled fermions [TVll4| preserve an exact chiral symmetry for a degenerate doublet 
of quarks, thereby realizing the minimal doubling of fermion species allowed by the Nielsen- 
Ninomiya theorem [15H17j . At the same time they also remain strictly local, and may thus be 
regarded as a cost-effective realization of chiral symmetry at non-zero lattice spacing, which 
is particularly suited for simulating a degenerate light doublet of up and down quarks. 

In a previous article [llj we began an investigation into the renormalization properties 
of a particular realization of minimally doubled fermions, described by the Borigi-Creutz 
action [3|Hll9llT0] . based on perturbation theory at one loop. In this paper|^ we present a similar 
analysis for another member of this class of fermions, proposed a long time ago by Wilczek [2], 
following previous work of Karsten [T]. Here we take the opportunity to revise and sharpen 
some of the conclusions of our earlier article pU] , by interpreting our results in a more general 
field-theoretical framework. In particular, what we interpreted in ref. as a shift in the 
quark's four-momentum can, in fact, be absorbed by adding appropriate counterterms to the 
lattice actions, so that any momentum shift disappears in the fully renormalized theory. The 
main focus of the present article is on the one- loop structure of minimally doubled fermions, 
taking the Borigi-Creutz and Karsten- Wilczek actions as particular examples. We identify 
the possible counterterms which can arise and study the consequences that can be inferred 
from their presence. 

Our findings can be summarized by the following: a consistent renormalized theory for 
minimally doubled fermions can be constructed by fixing the coefficients of three counterterms 
in the action which are allowed by the symmetries. In this work we determine these coefficients 
in perturbation theory, and discuss possible renormalization conditions to fix their values non- 
perturbatively. One of the principal results, namely the full expressions for the renormalized 
actions including all counterterms, can be found in eqs. (I75p and (I76p . We also present the 
computation of the vacuum polarization for both particular formulations. Our calculation 
demonstrates that radiative corrections at one loop do not introduce new divergences for this 
quantity. 

This article is organized as follows. After defining the two actions and the corresponding 
propagators and vertices in Section [2l we introduce in Section [3] the necessary counterterms 
which render the theories consistent under renormalization. With these tools in hand we fix 
the coefficients of the counterterms which appear in the quark action in SectionlU compute 
the matrix elements of quark bilinears and derive the conserved currents in Section[5] and 
also show that their renormalization constant is one. In Section [6] we present the calculation 
of the vacuum polarization of the gluon, while in Section [7] we discuss the implications of 
our findings for the practical implementation of minimally doubled fermions in numerical 
simulations. Finally, we make some concluding remarks. 



2 Actions, propagators and vertices 



In order to make this article self-contained, we recall here the basic definitions for Borigi- 
Creutz fermions [HHUHKlO] and also introduce the Karsten- Wilczek action [HIS]- In position 
space the lattice action of Borigi-Creutz fermions reads 



r 1 

X L fi=l 



^Preliminary results were reported in [12] , 
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- 'ip{x + afi) (7^ - [/T(x) '4}{x) + il){x) [mo - il;{x) 

where the matrices V and 7^ are defined by 

1 4 

r=2E>> 7; = r7Mr = r-7^, (2) 



with = 1. The Karsten-Wilczek action readsH 

4 

2a 



r 1 

'S'kw = ^ [V'Ca:^) (7/^ - ^74 (1 - <5m4)) + «/^) (3) 

'4){x + a/x) (7^ + 174 (1 - 5^4)) ii{x) + (^mo + ^/'(x 



We remind the reader that one Dirac spinor, ijj{x), in these expressions describes a degenerate 
doublet of quarks. In momentum space the free Dirac operator of Borigi-Creutz fermions 
takes the form 

2ir 

Vbc{p) = D{p) + D{p) + mo, (4) 

a 

with D{p) and D{p) given by 

^4 ^4 

^ip) = - Z (7m sin ap^l), D{p) = - Z (7m ap/,). (5) 
^J.=l M=i 

This action has two doublers, corresponding to the two zeros (Fermi points) at api = (0, 0, 0, 0) 
and ap2 = (vr/2, 7r/2, 7r/2, 7r/2). For Karsten-Wilczek fermions the free Dirac operator reads 

.4 .3 -3 

T^Kwip) = - ^lfj.smapf^ + --f4^{l-cosapk)+mo = D{p) + --f4^{l-cos apk)+mo. (6) 

In the latter case the term proportional to 74, which is chirally invariant since it anticommutes 
with 75, removes 14 of the doublers of the naive fermion action D{p), and only the doubler 
whose pole lies entirely in the temporal direction survives. The Dirac operator 2?kw(p) 
of the Karsten-Wilczek action exhibits then only two Fermi points, which are located at 
api = (0, 0, 0, 0) and ap2 = (0, 0, 0, tt) and describe - like for Borigi-Creutz fermions - two 
degenerate fermion species of opposite chirality. 

In principle, the term proportional to 74 can be multiplied by some coefficient. A, without 
spoiling the chiral symmetry of Karsten-Wilczek fermions, or modifying the minimal number 
of doublers. However, it is conceivable that for general values of A the transfer matrix does 
not exist or that the theory presents some other kind of problem. In this work we stick to the 
case A = 1, which is the common choice in the literature. El 

The two actions described here and investigated at length in this article represent two 
particular realizations of minimally doubled fermions, which respect chiral symmetry at any 
finite lattice spacing, but are no longer symmetric under the full hyper-cubic group. The 



"^In this article we will denote 4-vectors with Greek indices and spatial 3-vectors with Latin indices. The 
temporal component is /x = 4. 

^This coefficient A is in many ways similar to the Wilson parameter of Wilson fermions, which is commonly 
set to one in simulations. 
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Borigi-Creutz action is only symmetric with respect to a subgroup which preserves the hyper- 
cubic positive major diagonal. The term proportional to 74 in the Dirac operator of Karsten- 
Wilczek fermions selects, instead, a different particular direction in Euclidean space, the 
temporal axis. 

The breaking of hyper-cubic symmetry allows for mixing with operators of various di- 
mensions. In particular, mixing with lower-dimensional operators may occur, which implies 
the appearance of power-divergent coefficients proportional to 1/a". In the following we will 
show that for these actions the breaking of hyper-cubic symmetry indeed generates linearly 
divergent counterterms, whose coefficients need to be determined via some physical condition. 

We will see in the remainder of this article that in several instances the same reasoning 
applies to both actions considered. However, we would like to stress here that, despite many 
similarities, these two realizations of minimally doubled fermions are not equivalent. The ways 
in which the two actions are constructed, and the mechanisms for the removal of doublers, 
are indeed qualitatively different. El 

Although the distance between the two Fermi points pi and p2 is the same in each 
of the two actions considered (i.e., P2 — Pi = t^'^/'J^), these two actions cannot in fact be 
transformed into each other by a 4-dimensional rotation. On the other hand, in this article 
we will show that many results turn out to be qualitatively similar, owing to the common 
feature of the breaking of the hyper-cubic symmetry by a fixed direction in four-dimensional 
space (determined by the two Fermi points). 

It is easy to verify that both actions satisfy 75-hermiticity, with all the advantages and 
simplifications that this property implies, especially in numerical simulations. Another inter- 
esting observation is that the Dirac operators can be written as 



^|E7/.(Vm + V;) -ia74EVfcvA 

I u=l k=l ) 



(7) 



^KW = E7m(V^ + V;) -m74EVfcVfc^, (8) 

where V ^il^{x) = ^ \U^{x) ip {x + ajl) —\Ij{x)] is the nearest-neighbour forward covariant deriva- 
tive, and V* the corresponding backward one. Thus, it becomes apparent that the two real- 
izations of minimally doubled fermions bear a close formal resemblance to Wilson fermions, 

^4ii.on = ^ I E 7m(V^ + V* ) - ar ^ V; | . (9) 
I /x=i (1=1 ) 

Moreover, this demonstrates the presence of dimension-five operators in all three cases. For 
Wilson fermions the dimension-five operator breaks chiral symmetry, while for minimally 
doubled fermions chiral symmetry is preserved at the expense of introducing operators which 
break hyper-cubic symmetry. 

For the derivation of the propagators and vertices of Borigi-Creutz fermions we refer 



^There is some freedom in constructing actions of the Borigi-Creutz or Karsten-Wilczek type. For instance, 
in the latter case one can replace (1 — cos apfe) with other trigonometric functions without spoiling the basic 
properties of this kind of action, especially the fact that there are two Fermi points exactly located at api = 
(0, 0, 0, 0) and ap2 = (0, 0, 0, tt). 
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to Here we only remind that the Borigi-Creutz quark propagator can be written as 

~^ [^^^ '^"^ "'P^^ ~ ^ ^^"^^ + aniQ 

S{p) = a J , (10 

sin^ apfj^/2 + sin ap^f sin^ — 5 ^ sin^ app/2\ + {arriof 

the quark-quark-gluon vertex is given by 

I N ■ ( a{Pi+P2)ii , . a{pi+p2)t,\ , 
Vi{pi,P2) = -m ( 7m cos ^ -7/.sin ^ 1 (11 



(where pi and p2 are the incoming and outgoing quark momenta at the vertex), and the 
quark-quark-gluon-gluon vertex is 



V2ipi,P2) = ^iago [^7^ sin + 7^ cos j . (12) 

The fermion propagator of Karsten-Wilczek fermions is obtained by inverting the Dirac 
operator of eq. ([6]), and is given by 

4 3 
~^ sin apf, - 2i 74 ^ sin^ apk/2 + amo 

n / \ fi=l k=l 

'^^P^ = " 3 71 wl X • 

sin^ ap^ + 4 sin ap4 ^ sin^ o.pk/2 + 4 ^ sin^ apk/2 ^ sin^ a,pi/2 + (amo)^ 
At=l fc=l \fc=i / \i=i ) 

(13) 

As in the case of Borigi-Creutz fermions, but unlike many standard fermionic discretizations, 
we find that the denominator of this propagator does not possess a simple behavior under 
all momentum inversions (in this case of the fourth direction, so that it amounts to a time 
reversal) . 

After making the substitution vr+ap4, one obtains the propagator for the fermionic 

mode associated with the other Fermi point, the one at ap = (0, 0, 0, vr): 

4 3 

i ^ 7^ sin apfi_ — 2i 74 ^ sin^ aPfc/2 + arriQ 

'^yP) - ^ ^ 3 /3 \/3 \ ' 

sin^ ap^ — 4 sin ap^ ^ sin^ apk/2 + 4 ^ sin^ o.pk/2 ^ sin^ api/2 + (amo)^ 
/i=l k=l \k=i ) \i=\ ) 

(14) 

where we have also introduced a new set of Dirac matrices, 7^ = — 7^ and 74 = 74. If we now 
invert the direction of the four-momentum in this expression, the propagator of eq. ()13p 
is recovered. Since 75 = —75, we conclude that the modes corresponding to the two Fermi 
points have, as expected, opposite chirality. Note that the symmetry ap^ ^ vr — ap4 exchanges 
the zeros. 

In analogy to Borigi-Creutz fermions, we have also derived the quark-quark-gluon vertex 
for Karsten-Wilczek fermions, which reads 

I . ■ ( a{pi+p2)f, n r X • a{Pi+P2)A nK\ 
Vi{pi,P2) = -igo (7m cos ^ + 74 (1 - -5^4) sm ^ j , (15) 

while the quark-quark-gluon-gluon vertex comes out as 

T// ^ ^- 2f ■ «(pi+P2)m e X a{pi+p2)i,\ , . 

y2{pi,P2) = 2^«5o ( l/xsm - 74 (1 - V) cos I . (16) 
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The expressions for the vertices can also be easily derived by comparing the Dirac operator 
of Karsten and Wilczek with the Wilson-Dirac operator 



Dw{p) = - Ujf,smap^ + r {1 -cosap^)| +mo, 



(17) 



and observing that the hopping terms of these two actions are related by the replacement 
r o ^74, with a further restriction of the term proportional to 74 to its spatial components 
only. 

3 Counterterms 

A key result of this paper is the observation that the actions discussed in the previous section 
do not contain all possible allowed operators that are invariant under the subgroup of the 
hyper-cubic group which is left as a symmetry (preserving the positive major diagonal for the 
Borigi-Creutz case or the temporal axis for Karsten- Wilczek) . We will describe in detail how 
radiative corrections generate new contributions whose form is not matched by any of the 
terms in the original bare actions. These operators must then be added as counterterms in 
order to construct a consistent theory under renormalization, and this consistency requirement 
will uniquely determine the magnitude of their coefficients. 

In the following we will consider the massless case, ttiq = 0. Our task is to construct 
and add to the bare actions presented in Section[2] all possible counterterms which are allowed 
by the remaining symmetries. For this purpose we consider operators of dimension four or 
lower. We first classify them using the common notation in the continuum and then proceed 
to specify convenient lattice representations of these operators. 

We begin with the fermionic part of the actions. The presence of a conventional chiral 
symmetry strongly restricts the number of possible counterterms. Since they have to anti- 
commute with 75, we can restrict the list to those operators which contain 7^. Other Dirac 
matrices like 1, 75, 7^75 and a^^ can be excluded. The particular way in which symmetry 
breaking occurs in Borigi-Creutz fermions implies that we are allowed to construct operators 
where summations over single indices occur, in addition to the standard Einstein summation 
over two indices. Then also operators containing 7/^ = T are permitted. 

As a consequence, for Borigi-Creutz fermions there is only one possible counterterm of 
dimension four, namely tp T J2fi D^ip (which amounts to a renormalization of the speed of light 
for the fermions, relative to the positive diagonal axis). We can represent it on the lattice by 
writing it in a form similar to that of the hopping terms already present in the action. More 
precisely, we use the gauge invariant expression 



which is already present in the bare Borigi-Creutz action, albeit with fixed coefficient —2/ a. 
In the general renormalized action the coefficient of this operator must be tuned. It is con- 
venient in our perturbative work to use the convention that the operator ^lJ{x)T'^p{x) in the 
renormalized action has as a coefficient (—2 + C3)i/a, although for Monte Carlo simulations 



£4(50) Y XI (V'(^) r Uf,{x) ip{x + a/i) - %Ij{x + a/i) T C/^(x) V(x)) . 



(18) 



There is also one counterterm of dimension three, which is constructed from F, i.e. 




(19) 



a 



5 



other choices might be more appropriate. Since this piece of the action determines the lo- 
cation of the poles of the propagators, which are moved by radiative corrections, a possible 
renormalization condition is the requirement that the value of the coefficient must restore the 
poles to their original positions at pi and p2- 

For Karsten-Wilczek fermions things work out in a similar way. Here we are allowed 
to construct objects in which Kronecker deltas can constrain any Lorentz index to be equal 
to 4. It is easy to see that the only gauge-invariant counterterm of dimension four that can 
be added to the bare action is ip^iDi^ip. A suitable discretization for this operator is 

The counterterm of dimension three, i.e. 

^^3(50) 



diigo) Y (V'(^) ^4 Ui(x) i^{x + aZ) - tPix + aZ) 74 ul{x) V(x)) . (20) 



■V(x) 74^^(2;), (21) 

is already contained in the bare Karsten-Wilczek action, where it has a fixed coefficient of 3/a. 
In the fully renormalized theory we will denote the coefficient of the term ■i/'(x)74'(/'(x) by 
(3 + (is) i/a. 

In perturbation theory the coefficients multiplying the above counterterms are functions 
of the coupling, starting at order Qq. They generate new vertices and propagator insertions. 
At one loop they give rise to additional contributions to fermion lines, and these insertions 
must be taken into account for a consistent one-loop calculation. From the above expressions 
the rules for the corrections to external fermion propagators, which will be needed for the 
calculations presented in this paper, can be easily derived. These are independent of the 
lattice discretization chosen for the counterterms. Since the propagator is the inverse of the 
quadratic part of the action, they are given in momentum space by 



for Borigi-Creutz fermions, and 



^cMrJ:p., (22) 



■r / N ^^^3(50) ,^o^ 

««4(5o) 74P4, 74 (23) 

a 



for Karsten-Wilczek fermions, respectively. We will determine all their coefficients at one loop 
in perturbation theory (see Section Hj) , by requiring that the renormalized self-energy assumes 
its standard Lorentz-invariant form. 

We also need counterterms for the pure gauge part of the actions of minimally doubled 
fermions. Although at the bare level the breaking of hyper-cubic symmetry is generated by 
the fermionic actions only, it propagates via the interactions between quarks and gluons also 
to the gauge sector in the renormalized theories. One effect is that some of the terms in the 
purely gluonic part can renormalize with different factors, and as a consequence pure gauge 
counterterms must be added to the renormalized actions to correct this imbalance. They are 
of the (continuum) Tr FF form, but with non-conventional choices of the indices which reffect 
the breaking of the Lorentz (hyper-cubic) symmetry. 

Let us consider first the Borigi-Creutz case. If we choose all four indices of TrFF to 
appear only once in each summation (as allowed by hyper-cubic symmetry breaking), we can 
construct a counterterm which has the continuum form 

^ TvFxp{x)F^r{x). (24) 

X,p,(7,T 
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However, since F^^, is antisymmetric in the indices, this is identically zero. The next possi- 
bility is to contract one Lorentz index shared by both field tensors (as in the usual Einstein 
convention), whilst summing individually over the two remaining indices, i.e. 

CP (go) E TrFAp(x)Fp,(x). (25) 

A,p,r 

This operator, whose form is reminiscent of the energy-momentum tensor, is the only possible 
purely gluonic counterterm for this action. H A lattice counterpart for this counterterm can 
be obtained by employing the widely used "clover" expression of the F^,^ tensor [18j. 

At one loop this counterterm contributes only via insertions in gluon propagators. De- 
noting the fixed external indices at both ends of these lines with fi and i^, all possible lattice 
discretizations of this counterterm yield the same Feynman rule in momentum space, namelyH 



cp(ffo) 



A A 



(28) 



As we will see explicitly in Section [6l the presence of this counterterm is essential in order to 
ensure the correct renormalization of the vacuum polarization. 

It is not hard to infer that in the case of Karsten-Wilczek fermions the temporal plaque- 
ttes (the chromo-electric field) renormalize differently compared with the spatial plaquettes 
(corresponding to the chromo-magnetic field). The counterterm to be introduced will contain 
an asymmetry between these two kinds of plaquettes, and can be written in continuum form 
as 

dpigo)Y.TrFpxix)F^xix)Sp^. (29) 

p,A 

This is the only required gluonic counterterm for this action, since introducing another factor 
of 6\4 in the above expression will produce a vanishing object. It is immediate to write down 
a lattice representation for it, using the plaquette: 

^p(ffo)f E(i-]^^^4a(x)). (30) 

p,A 

The Feynman rule for an external gluon line that this counterterm generates reads 

- dp{go) p^py (5^4 + 6yi) - 6fj,„ [p^ S^iSyi + pf) ■ (31) 

This term guarantees the correct renormalization of the vacuum polarization. 

From the gauge-invariant expressions of the operators introduced in this section we can 
read off that interaction vertices are generated by the counterterms. However, these vertex 



^In fact, choosing two pairs of summed indices will reproduce the usual Lorentz invariant term 

^Tt F\p{x) F\p{x), already present in the bare action. 
®The lowest order in the expansion of Fxp(x) Fpt(x) in momentum space is 

- pxPpApAr - PrPpApAx + plAxAr + PxPrApAp. (26) 

When the external indices of the gluon propagator are set equal to ^ and v, this becomes 

- P^iPxAp,A^ - Pt^pxApAy + p^ApAt, + pxp^-ApAy Sp^,. (27) 

Here all indices but ^ and v are summed. A similar derivation holds for the Karsten-Wilczek action, where it 
gives eq. ((3T|) . 
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insertions are of higher order in (at least of 0((7g)), and thus they cannot contribute to the 
one-loop amplitudes presented in this paper. 

It is worth stressing that the form of the counterterms that we have constructed remains 
the same at all orders of perturbation theory. Only the values of the coefficients change 
according to the loop order considered. The same counterterms appear also at the non- 
perturbative level, and will be required for a consistent numerical simulation of these fermions, 
as we will see in Section [71 Finally, we want to emphasize here that counterterms not only 
provide additional Feynman rules for the calculation of loop amplitudes, but can also modify 
Ward identities and hence, in particular, contribute additional terms to the conserved currents, 
as we will see explicitly in Section [5l 



4 Fermionic counterterms at one loop 

We now show how the coefficients of the counterterms that appear in the quark action can 
be fixed by computing the quark self-energy. Technical details in the case of Karsten-Wilczek 
fermions are deferred to appendix lA.ll while the case of Borigi-Creutz action was described 
at length in pT] . 

In the following we will work in some general covariant gauge, where d^Af^ = and a 
denotes the gauge parameter. For Borigi-Creutz fermions, the result for the one-loop diagrams 
of the quark self-energy, without including the counterterms, reads (see eq. (15) of [llj) 

S(p,mo) = i/Si(p) +moS2(p) + ci(fifo) • iF Vp^ + 02(50) • i -, (32) 



where 



log aV + 6.80663 + (1 - q) ( - log aV + 4.79201o) 
4 logaV - 29.48729 + (1 - q)( - logaV + 5.79201o) 



(33) 
(34) 



ci(5o) = 1.52766 -^Cf, (35) 
C2(5o) = 29.54170 -^Cf, (36) 

with Cp = {N^ — l)/2Nc. From eq. ([32]) alone, the inverse propagator at one loop level would 
assume the form 

S-i(p,mo) = (l-Si) .{i^+mo(l-S2 + Si) -ici(5o)r5]p^-^^^r}. (37) 

We now show how the coefficients multiplying the counterterms of the Borigi-Creutz quark 
action can be fixed by imposing the condition that the renormalized propagator take the 
standard form 

Z2 

S(p,mo) = , (38) 

tp + Am mo 

where the wave-function and quark mass renormalization factors are given by 

Z2 = (l-Si)"\ (39) 

Zra = 1-(S2-Si), (40) 
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Indeed, by looking at eq. ([2^ we see that the piece in eq. ([57]) which is proportional to 
ci{go) can be eliminated by tuning the coefficient of the counterterm of dimension four, 
04(50) tpT J2fj. D^tp, while the one proportional to 02(50) can be eliminated by a suitable choice 
of the counterterm of dimension three, ic^^go) ipTip/a. After inserting these counterterms, one 
has 



ip + mo 



+ 



1 



+ mo 



Si + mo S2 + ci • i r ^ + C2 • i Ci-iT^Pfj, 

1 



r 

C3 • z - 

a 



if + mo 



1 



'(1 - Si) + mo (1 - S2) - (ci - Ci) ■ iVY^^p^- (c2 - C3) • i - 



r • 



(41) 



It follows that 04(50) = ci(5o) and 03(50) = £2(50)1 if eqs. ([55j) - ([in|) are to be recovered. Thus, 
at this order the coefficients of these counterterms for Borigi-Creutz fermions are determined 



as 



£3(90) 
£4(90) 



29.54170 



1.52766 



167r2 
9l 



Cf + 0(5, 



167r" 



Cf + 0(5o'). 



(42) 
(43) 



Of course, in order to be useful for numerical simulations, these coefficients will have to 
be determined non-perturbatively, by imposing suitable renormalization conditions. We will 
return to this question in Section [71 

The explicit calculation of the self-energy for Karsten-Wilczek fermions at one loop pro- 
ceeds along the same lines as for Borigi-Creutz. Further details as well as the individual results 
for the tadpole and sunset diagrams can be found in Appendix lA.il The result for the total 
contribution, without counterterms, is 



S(^),mo) = i^Si(p) + mo'E2{p) + (ii(5o) -i^APi + ^2(50) ■ i — , 



where 



S2(P) 

diigo) 
^2(50) 



167r2 
167r2 



Cf 



log a V + 9.24089 + (1 - q) ( - log aV + 4.792010 
4 logaV - 24.36875 + (1 - q)( - logaV + 5.79201o) 



-0.12554 • 



-29.53230 



167r2 
■A 



Cf, 



Cf. 



Without any counterterms, the inverse propagator at one loop is 

S-i(p,mo) = (l-Si) • (i/ + mo(l-S2 + Si) - idij4P4 



■74 



(44) 

(45) 
(46) 

(47) 
(48) 

(49) 



As in the case of Borigi-Creutz fermions, the extra contributions proportional to di(5o) and 
^2(50) can be cancelled by suitably tuning the coefficients of the counterterms ^4(50) ip^fiD^tp 
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and id^{gQ) ip^^tp/a respectively. Thus, the full inverse propagator at one loop for Karsten- 
Wilczek fermions can be written in the standard form, eqs. ()38p - ()40p . The coefficients of the 
counterterms for Karsten-Wilczek fermions so determined are at one-loop order 

^3(50) = -29.53230 •^C7f + 0(5^), (50) 

diigo) = -0.12554 •^Cf + 0(5o'). (51) 

One may expect that the above subtraction procedure can be carried out consistently at every 
order of perturbation theory. After the subtractions via the appropriate counterterms are 
properly taken into account, the extra terms appearing in the self-energy can be eliminated. 

Recalling that the term proportional to 74 in the Karsten-Wilczek action can be multiplied 
by a parameter A, one may wonder whether a suitable choice of A could eliminate the power- 
divergent contribution to the quark self-energy, without resorting to counterterms. This 
could in principle be accomplished in view of the fact that the tadpole and sunset diagrams 
contribute with opposite sign to the self-energy, while the tadpole contribution is linearly 
proportional to A. We have investigated this issue using our perturbative expressions, and 
concluded, however, that such a cancellation cannot take place. The reason is that while 
the value of the tadpole decreases as A is lowered, the contribution of the sunset diagram 
also decreases at the same time. Since the latter diagram always remains much smaller than 
the tadpole it cannot compensate its value. Thus, it is not possible to eliminate this power- 
divergent extra term without counterterms. 



5 Quark bilinears and conserved currents 

For many applications, knowledge of the renormalization factors for quark bilinears is required. 
These are obtained by computing the appropriate vertex diagrams and adding the contribution 
Si to the quark self-energy. For Borigi-Creutz fermions, the results for the vertex diagrams 
of the various bilinears are given in [llj. For Karsten-Wilczek fermions the corresponding 
results are listed in Appendix lA.2[ Here we note that for the latter, hyper-cubic symmetry 
breaking induces different radiative corrections for temporal and spatial components of the 
vector and axial- vector currents. 

In the following we discuss the issue of the renormalization of the quark mass. For both 
realizations of minimally doubled fermions, chiral symmetry is preserved and so the quark 
mass does not undergo any additive renormalization. Thus, the relation between the bare 
and renormalized quark masses is 

mR = Zmmo, (52) 

where is given in eq. ()40p . The full expression for the renormalization factors of the scalar 
and pseudo-scalar densities in perturbation theory at one loop is 

Zs = Zp = 1-(As + Si), (53) 

where As is the vertex correction of the scalar density, which is given by eq. (I82p for Karsten- 
Wilczek and eq. (25) of pT] for Borigi-Creutz fermions, respectively. This number is exactly 
equal to the 0((7g)-contribution to S2 (eq. ([Mj) or ([iG]) !. but comes with an opposite sign: 
As = -5^2- Thus, when we compare them with eq. (j40p . we see that the renormalization 
factors Zs and Zp of the scalar and pseudo-scalar densities satisfy 

1/Z^ = Zs = Zp, (54) 
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where the last equahty is a consequence of chiral symmetry. We have thus verified at one loop 
that the renormalization of the quark mass for both variants of minimally doubled fermions 
presented here has the same form as, say, in the case of overlap or staggered fermions. 

Using the expressions for the vertex corrections of the local vector and axial-vector cur- 
rents and taking the renormalization of the wave- function into account, the corresponding 
renormalization factors and Za are not equal to one. In order to construct the conserved 
currents, which are protected against renormalization, one has to derive the chiral Ward iden- 
tities, for instance, along the lines of ref. [19]. To do this, we have to use the expressions for 
the lattice actions of Borigi-Creutz and Karsten-Wilczek fermions in position space, which are 
given in eqs. ([TJ and ([3]) respectively. 

It is important to note that the counterterms which have been added to the actions can 
also contribute to the chiral Ward identities, and thus provide new terms in the expressions 
of the conserved currents. Actually, it is easy to see that the counterterms of dimension three 
do not modify the Ward identities, but the counterterms of dimension four do. The latter 
generate additional terms in the Ward identities and hence also in the conserved currents. 
Thus, the form of the conserved currents turns out to be different from what the bare actions 
would have given. 

It is worth stressing again that, as we already discussed in |11) . these actions in the 
massless case are invariant under a chiral U(l) U(l) transformation. The chiral Ward 
identities associated with these exact symmetries yield the currents of the theory. The chiral 
U(l) ® U(l) symmetry of minimally doubled fermions implies that as the quark mass is tuned 
to zero there is only one Goldstone boson, which can be naturally considered to be the neutral 
pion. The charged pions will instead be massive in the chiral limit (at non-zero lattice spacing). 

If one applies the standard vector and axial transformations, i.e. 



under which the Lagrangian remains invariant, one can identify the conserved vector current 
for Borigi-Creutz fermions in the renormalized theory as 



6vip = iayip, 



Svi' = —iipay, 
5 Alp = iV'ttATs, 



(55) 




(56) 



The axial-vector current (which is only conserved in the massless case) is given by 




i^ix) (7^ + i 7|,)75 Uf,{x) ipix + ap,) + ipix + ap,) (7^ - i 7|,)75 ?7^(x) ipix) 
(lp{x) r75 U^{x) iP{x + ap) + V^(x + a/i) r75 UUx) ^P{x)] . 



(57) 



For Karsten-Wilczek fermions the conserved vector current turns out to be 
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and the axial-vector current is given by 

Al,{x) = ^ (^i;{x) (7^ - i74 (1 - (5^4))75 U^{x) i){x + a/I) 

+tl;{x + afi) (7^ + i74 (1 - (5^4))75 U^ix) ^^{x) 

(^^(^) ^4^5 Ui{x) il,{x + a4) + + a4) 7475 f/i(x) i^{x)^ . (59) 

It is instructive to verify that the renormahzation constants of these currents is equal to one, 
since this exercise illustrates the important role of the counterterms. Here we restrict a detailed 
discussion to the conserved vector current for Karsten-Wilczek fermions, noting that this is 
entirely analogous to the Borigi-Creutz case. Furthermore, the corresponding expressions 
for the conserved axial-vector current in both discretizations are trivially obtained from the 
formulae below by replacing 7^ by 7/^75 and 74 with 7475. 

The renormahzation factor of the vector current is given by 

Zv = 1- (Av + Si), (60) 

where Ay denotes the vertex correction and Si the self-energy. The sum of the vertex (dia- 
gram (a) in Fig.[l]), the "sails" (diagrams (b) and (c)) and the operator tadpole (diagram (d)) 
corresponds to the first two lines in eq. (f58l) and yields 



^0 r' o, 



logaV - 9.24089 + ■ 0.12554 + {I - a)(logaV - 4.79201 



(61) 



Now, a counterterm is also included in the expression of the conserved current (c.f. the last 
line of eq. (j58|) ). and contributes a factor of 

^4(50) 74 = ^4(50) 7m'^m4 (62) 

to its renormahzation factor. At lowest order in perturbation theory this was evaluated as 
part of the determination of the self-energy, with the result 

d4(5o) = -0.12554 -^^0^ + 0{gt). (63) 

This cancels exactly the contribution proportional to 5^4 in eq. ()6ip . So, apart from including 
the wave-function renormahzation, the result for the proper diagrams is 



9l 



lOTT^ 



logaV - 9.24089 + {I - a)(logaV - 4.79201 



(64) 



Finally, we see that this expression is equal and opposite to the contribution of Si(p) of 
the quark self-energy, eq. ()45p . and so exactly compensates it. Thus, we conclude that the 
renormahzation factor of these point-split currents is one. This holds to all significant digits 
that we have achieved, and confirms that the expressions that we have derived via the chiral 
Ward identities are indeed conserved currents. 

For Borigi-Creutz fermions, the sum of vertex, sails and operator tadpole diagrams for 
the conserved vector current is pij 



^2 



- log a V - 6.80664 (1 - a) ( log a V - 4.7920l) 



.,2 



1.52766 --^Cp-r. (65) 
167r^ 
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Again, at this point one has to add the contribution from the counterterm of dimension four 
in the conserved current (i.e. the last line of eq. ([56]) ). 

C4(5'o)r, (66) 

whose coefficient was already fixed by the result of the one-loop self-energy: 

C4(ffo) = 1.52766 --§^0^ + 0(9o)- (67) 

This cancels the Lorentz non-invariant term in eq. (|65p . After including the contribution of 
the wave-function renormalization it can be easily seen that the renormalization constant of 
this current is equal to one as well. 



6 Vacuum polarization 

In order to gain a deeper understanding of the properties of minimally doubled fermions, we 
have also taken the task to calculate the one-loop vacuum polarization of the gluon for our 
two realizations of minimally doubled fermions. 

Here we focus on the radiative corrections to the bare gluon propagator which arise from 
fermion loops and the gluonic counterterm. At one loop the perturbative contributions to the 
vacuum polarization due to loops of gluons and ghosts are independent of the chosen fermionic 
lattice action, and are thus irrelevant for the problem we are studying. However, quark loops 
are able to generate hyper-cubic-breaking terms, and here we show that this is what indeed 
happens for both Karsten-Wilczek and Borigi-Creutz fermions. 

It is instructive to recall the expression for one flavour of Wilson fermions. In this case, 
neither breaking of hyper-cubic symmetry nor fermion doubling take place, and the result is 



167r2 



Co 



- log p^a^ + 4.337002 



(68) 



where C2 is defined via Tr {tH^) = C2 S"^. It is straightforward to see that this gauge invariant 

( f) 

result satisfies the Ward identity p'^HfiJ {p) = 0, which expresses the conservation of the 
fermionic current. 

The result of our calculations for the above quantity using Borigi-Creutz fermions, without 
including the gluonic counterterm, reads 



PfiPiy - ^t^iyP 



^C2(-^logpV + 23.6793 



(Pm + p.) E pa - / - 'J/.. ( E pa) ' ) t|^C2 • 0.9094 



(69) 



By comparing the coefficient of the divergence to the expression for Wilson fermions we 
see explicitly that this result corresponds to two flavours of quarks: Each of the two doublers 
contributes an equal amount to the divergence. The same is true for Karsten-Wilczek fermions, 
for which our calculation gives the result (again without counterterms) 



ni/J(p) 



PfiPu - S^uP"^ 



167r2 



■logp^a^ + 19.99468 



if b^^b,^ + p1) ) -§^C2 ■ 12.69766 



(70) 



167r2 
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We first note, looking at the second lines of the above equations, that additional terms appear, 
compared with a standard situation like Wilson fermions. A remarkable observation is the 
fact that although each of the two actions breaks hyper-cubic symmetry, the extra terms 
still satisfy the Ward identity p^HfiJ (p) = 0, expressing current conservation, as can be easily 
verified. It is easy to see that one cannot obtain smaller expressions for the symmetry-breaking 
terms which are symmetric in fj, and and still conserve the current. Moreover, our one loop 
results are the most general symmetric quadratic functions of p which can be constructed 
using the hyper-cubic breaking objects 6^4 or J2xPx- 

The correct renormalization of the polarization of the vacuum requires the inclusion of 
the counterterm of the pure gauge action. To see how this works in detail, we consider for 
illustration Borigi-Creutz fermions. Adding the one-loop results to the tree-level expression, 
and exploiting pf^HfiJ (p) = in the usual way, we can write 



Oitu — — — — Ot T — r „ ^ 



p. p. 



T T 



Spu — a 



p'^ 



p2 



= pHi-uih) ^ ^ ^''^ 

The function n(p^) in the first line of the above equation multiplies the standard Lorentz 
invariant one-loop expression, which is then rearranged such as to produce the functional 
form of the continuum tree-level gluon propagator in the last line. From this one can read 
off that the gauge parameter is renormalized according to 5 = a (1 — 11) + 11, and that 
Z3 = 1/(1 — 11(0)). The remaining terms, the ones proportional to n(p^), are those that break 
hyper-cubic symmetry, and they cannot be rearranged in a similar way. 

It is thus evident that these hyper-cubic-breaking contributions must be eliminated, and 
this can be achieved by employing the gluonic counterterms which we introduced in Section[3l 
Indeed, the expression for the gluonic counterterm for Borigi-Creutz fermions in momentum 
space, eq. ()28p . is structurally identical to the additional terms in the vacuum polarization. 
Requiring the one-loop vacuum polarization to assume the standard Lorentz invariant form 
(that is, only the first term in the above result) then uniquely determines the coefficient of 
the counterterm. The non-standard contributions are cancelled for Borigi-Creutz fermions if 
we set ^ 

cp(<7o) = -0.9094 --^02 + Oig^). (72) 
The reasoning for Karsten-Wilczek fermions is entirely analogous, and one concludes that 

dp{go) = -12.69766 ' C2 + 0{g^). (73) 

The most important thing to realize is that there are no power-divergences in our results for the 
vacuum polarization. In principle such divergences could arise with coefficients proportional to 
1/a^ or 1/a. We have explicitly checked in our calculations that the l/o^ tadpole contribution, 
when non-zero, is in all cases of equal magnitude and opposite sign with respect to the sunset 
diagram. 



^It is interesting to note that the numbers for these diagrams are much larger than in the case of Wilson 
fermions, where the coefficient of (;oC2/167r^ for the tadpole is —9.67590. For Karsten-Wilczek fermions this 
number turns out to be —36.31464 for each spatial component and 7.12931 for the temporal component. For 
Borigi-Creutz fermions it is even larger, —73.71980. 
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We can understand on general grounds why such power-divergences cannot appear. To 
construct hyper-cubic breaking terms one has to employ objects like F and J2fiPiJ, i^^^ Borigi- 
Creutz fermions) and 74 and p4 (for Karsten-Wilczek fermions). However, after the traces 
of the fermions loops are evaluated no Dirac structures are left over, and momenta cannot 
anyway appear at the 1/a^ level. Linear pieces in the momenta, which would be required in 
case of a 1/a power divergence, are instead prohibited by the symmetry of the diagrams. 

We have also discovered that the hyper-cubic-breaking terms can be put for both actions 
in the same algebraic form, namely 

pHi^,T}{^u,r} + 6^A^,T}y,T} - \ {^,r}({7^,/}{7,,r} + {i.Mi^.^}). (74) 

where in the case of Karsten-Wilczek fermions F must be replaced by 74/2. This substitution is 
suggested by comparison of the standard relation = \ J2 ij.il t^+lp) of Borigi-Creutz fermions 
with the formula 74 = ^ J^fiilfi'^l'^j.) Karsten-Wilczek fermions, expressing the symmetries 
of the action (as can be seen from Section [21 when one expands the propagator of the latter 
action around the second Fermi point). Whether there is any deeper significance to this 
structural "equivalence" of the hyper-cubic-breaking structures in the vacuum polarizations 
remains an open question. 



7 Numerical simulations 



In this section we discuss the implications of our one-loop perturbative calculations for nu- 
merical simulations of minimally doubled fermions. The first thing to note is that simulations 
must be based on the complete renormalized actions, including the counterterms. In position 
space the full expression for the Borigi-Creutz action reads 

Sbc = ^'^^\^Y1 k(^) + C4(5o) r + i7^) Uf,{x) tpix + ajl) 

- ijj{x + aJT) (7^ - Ci{go) F - i7^) [/^(x) ^(x) 



+ V(x) 7710 + 03(50) — ]ip{x) (75) 

+ /^E (1- l^I^eTrP^J +cp(5o) E '^F^.p{x)Fpu{^ 

where we have redefined the coefficient of the dimension-three counterterm, via C3(g'o) = 
—2 + 03(50)) and F is some lattice discretization of the field-strength tensor. The renormalized 
action for Karsten-Wilczek fermions reads 

r ^ 4 _ 

'S'kw = E ^ ^ E k(^) + ^4(50) 5pi) - i74 (1 - 5^4)) U^{x) ^(x + a/2) 
- %l){x + afi) (7^(1 - ^4(50) '^^4) + ilA (1 - S^a)) Ul{x) tpix) 
+ V^(x) (mo + ^3(50) ^) ^{x) (76) 
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where d^{go) = 3 + d3((7o)- perturbation theory the coefficients of the counterterms admit 
the expansions 

Mao) 
CP (go) 

and 

^3(50) 
diigo) 
dp (go) 

In order to define the fuh actions at the non-perturbative level one must impose suitable renor- 
malization conditions which fix the values of the counterterms beyond perturbation theory. 
Below we discuss possible scenarios how this could be achieved. 

It is fairly straightforward to determine the coefficients of the counterterms of dimension 
four via a non-perturbative condition: as we have seen in Section[5l these counterterms ensure 
that the conserved currents have unit renormalization. The coefficients C4 and ^4 can then 
be fixed by requiring that the electric charge be equal to one. To this end one can compute 
suitable ratios of three-point and two-point correlation functions, involving the expressions in 
eqs. (f56]l and (|58l) . respectively. Adjusting the coefficients until the charge is unity fixes the 
values of C4 and d^. It is an empirical fact that ratios of correlators are obtained with good 
statistical accuracy. 

Furthermore, radiative corrections induce a shift of the poles of the quark propagator 
away from their tree-level positions. Provided that their coefficients are appropriately tuned, 
the counterterms of dimension three ensure that the two Fermi points can be moved back to 
their original locations. It is important to realize that radiative corrections, when moving the 
poles, do not introduce a sign problem into the Monte Carlo generation of configurations. The 
gauge action remains real, and the eigenvalues of the Dirac operator are always in complex 
conjugate pairs, making the fermion determinant always non-negative. 

On the other hand, these shifts can introduce oscillations as a function of separation into 
some hadronic correlation functions. Such oscillations, familiar from the staggered formula- 
tion, come about since the underlying fermion field can create several different species, and 
these species occur in different regions of the Brillouin zone. It would be interesting to explore 
whether or not these oscillations could be cancelled by constructing hadronic operators spread 
over nearby neighbours. 

It is important to remember that because the two species are of opposite chirality, the 
naive 75 matrix is physically a fiavour non-singlet. The naive on-site pseudoscalar field V'TsV' 
can create only fiavour non-singlet pseudoscalar states. To create the fiavour-singlet pseu- 
doscalar meson, which gets its mass from the anomaly, one needs to combine fields on nearby 
sites with appropriate phases. 

The purely gluonic counterterm for Borigi-Creutz fermions introduces new operators into 
the renormalized action, in which chromo-electric and chromo- magnetic fields enter. Since the 
positive diagonal of space time hypercubes is selected as special, terms involving J2ij.,u,p ^i^pPup 
can enter. This gives rise to combinations like E ■ B, E1E2, B2B3 and similar. These cross 
terms can be removed by a diagonalization process, essentially a rotation redefining the time 
direction to be along the positive diagonal. However, then the coefficents of E'^ and B^ can 



-2 + 4%2+4^o^ + ... 

c?9l + c?\'o + ■■■ 



(77) 



3 + 41)^2^42)^4^... , 
d^^gl + df^g', + ..., 
d?gl + d^P^g'o + .... 



(78) 
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differ. Effectively, the speed of light for the gluons has been renormalized. The coefficient 
cp could then be fixed by tuning its value to restore the E, B symmetry. The above effects 
could turn out to be small, given that at tree level only the fermionic actions break hyper- 
cubic symmetry. It could also happen that other derived quantities are more sensitive to this 
coefficient and more suitable for its non-perturbative determination. In general, one can look 
for Ward identities and study their deviation from the standard Lorentz invariant form, as a 
function of cp. 

For Karsten-Wilczek fermions the purely gluonic counterterm induces an asymmetry be- 
tween the plaquettes containing a temporal link relative to those involving spatial links only. 
Fixing the coefficient of this counterterm, dp, could then be accomplished by computing a spa- 
tial plaquette or Wilson loop, and then equating its result to its counterpart with components 
in the time direction. 

Eventually Monte Carlo simulations will reveal the actual amount of symmetry breaking, 
which could turn out to be large or small, depending on the observable considered. One 
important such quantity is the mass splitting of the charged pions relative to the neutral pion. 
Furthermore, the relative magnitude of these symmetry-breaking effects could be substantially 
different for Borigi-Creutz and Karsten-Wilczek fermions. Which of these two realizations has 
the potential to become the preferred choice for numerical simulations, will largely depend on 
this issue. 

So far we have not touched on the important subject of identifying the leading lattice 
artefacts associated with minimally doubled fermions. As we observed in connection with 
eqs. ([7]) and ([8]), both realizations of minimally doubled fermions contain a dimension- five 
operator in the bare action, and this leads one to expect the leading lattice artefacts to be 
of order a.H Naively, one might assume that the preservation of chiral symmetry would 
automatically ensure 0(a) improvement, but here we can see that this is not always the case. 
In minimally doubled fermion actions there are in fact manifest 0(a) effects, which arise as a 
consequence of the breaking of hyper-cubic (and not chiral) symmetry. 

This naturally leads to a discussion of the subject of 0(a) improvement. Since these 
actions are not improved from the beginning, at least one further dimension-five operator will 
be needed in order to cancel 0(a) contributions in on-shell matrix elements. Of course all 
possible operators of dimension five that are consistent with the symmetries must be con- 
sidered. Since minimally doubled fermions respect chiral symmetry, this excludes operators 
like the "clover" term, which contains a^u and therefore does not anticommute with 75. We 
attempted a cursory analysis, but there appear to be quite a number of possible operators, 
even after some of them are reabsorbed into a rescaling of parameters or eliminated using the 
equations of motion. Just to give an example, among the possible operators for Borigi-Creutz 
fermions we find 4'^Yl,^,v -^i^-^f'^ ^'^d Yl,^i,v,\^tiy^\^iJ.i^- We emphasize that such additional 
dimension-five operators can occur not only in the quark sector, but also in the pure gauge 
part. In fact, when Lorentz invariance is broken, the statement that only operators with even 
dimension can appear in the pure gauge action is no longer true. As this discussion shows, 
the issue of 0(a) improvement can be quite intricate for this particular type of fermionic dis- 
cretizations. The task of classifying the minimal set of independent operators is an interesting 
problem but beyond the scope of this paper. 

We close this section with some speculations. First, it is not entirely clear that we need 
to move the Fermi points back to the free theory positions. As long as there are still two poles 

*This contradicts the findings of the (numerical) tree-level analysis of ref. 5 , according to which cutoff 
effects for pion masses and decay constants were of order a? . 



17 



in the propagator, we have minimal doubUng. Further investigations could reveal whether 
relaxing this condition simplifies the tuning of counterterms. Finally, if we are willing to live 
with an induced anisotropy, we only need to tune C3 and C4 roughly and then use physical 
observables to measure the remnant anisotropy. Then only cp needs to be adjusted so that 
the anisotropy for the gluons is the same as for the fermions. 

8 Conclusions 

Following on the first investigation of minimally doubled fermions for Borigi-Creutz fermions 
we have presented further perturbative study including another realization, namely 
Karsten-Wilczek fermions. At the same time, we have re-interpreted our analysis of the 
renormalization of both these actions in the context of a solid field-theoretical framework. 

Our investigations show that both Borigi-Creutz and Karsten-Wilczek fermions are de- 
scribed by a fully coherent quantum field theory. Since the complete set of operators allowed 
by the symmetries of these fermions includes several not present in the tree- level actions, 
counterterms must be introduced for a consistent renormalized theory. After adding the 
counterterms, only a small number of new mixings arises for the matrix elements of local bi- 
linears, none of which is power divergent. For the local vector and axial-vector currents, some 
finite mixing with hyper-cubic breaking operators occurs. For the scalar and pseudoscalar 
densities and the tensor operator, on the other hand, the structure of the mixings is as in the 
continuum. 

We have also constructed the conserved vector and axial-vector currents for both kinds 
of fermions. The vector current is isospin-singlet, representing the conservation of fermion 
number. The axial current, however, is a non-singlet because the doubled fermions have op- 
posite chirality. These currents have simple expressions which involve only nearest-neighbour 
points, and do not undergo any mixing. We have verified at one loop in perturbation theory 
that their renormalization constants are equal to one. One of the most attractive features 
of Borigi-Creutz and Karsten-Wilczek fermions is that they belong to the very few lattice 
discretizations that yield a simple and essentially ultralocal expression for a conserved axial- 
vector current. These features could turn out to be a key advantage in numerical simulations. 
Furthermore, we have also calculated the polarization of the vacuum for both actions. We 
have proven that, in spite of the breaking of hyper-cubic symmetry, no power divergences 
appear for the vacuum polarization of the gluon. 

In summary, we have constructed the renormalized theory up to and including 0((7o) for 
the two realizations of minimally doubled fermions considered. We have discussed perturbative 
and non-perturbative conditions for fixing the coefficients of the three counterterms required 
for both realizations. We have argued that under reasonable assumptions and following the 
determination of these counterterms, no special features of these two realizations of minimally 
doubled fermions should hinder their successful Monte Carlo simulation. 

Some questions merit further consideration: firstly, one should revisit the problem of 
formulating conditions which allow for a precise determination of the coefficients of the 
dimension-three counterterm at the non-perturbative level, and the same applies to the glu- 
onic counterterms. Secondly, attempts to improve convergence towards the continuum limit 
must take account of the inherent hyper-cubic symmetry breaking and the induced mixing 
with dimension-five operators. 
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A One-loop calculations for Karsten-Wilczek fermions 

In this appendix we report the individual expressions at one-loop for various quantities com- 
puted using the Karsten-Wilczek action. Figured] lists all diagrams which are needed for the 
perturbative calculations presented in this article. 



A.l Self-energy 

Using the expression for the vertex V2 (p, p) of eq. (|16p , the tadpole contribution to the self 
energy can be easily computed. In a general covariant gauge, where d^A^ = 0, its expression 
is 

4 • ^ (1 - ^(1 - «)) • iaglCp^Yl ^^^^Pt^ - E 74(1 + 0(g2)) 

M=i k=i 

= ^o'Cpf (l-i(l-a))(z^-^)+0(a), (79) 



where Zq is given by [20H22] 

and terms of 0(a) and higher can eventually be set to zero. The term proportional to i;^ 
is the same as for Wilson fermions, while the other term would imply a power-divergent 
mixing of order 1/a with the dimension-3 operator il^^^ip, unless there be a cancellation with 
an analogous term coming from the contribution of the sunset diagrams to the self-energy 
(contained e.g. in diagram (e) of Figure[T|). Here we show that no such compensation occurs. 

We have computed the sunset diagram using special computer codes written in FORM 
[23l[24l and Mathematica, and also checked it against calculations by hand. The result is 



S™*(p,mo)=i^-^CF 



log a V - 2.99216 (1 - a) ( - log a V + 7.850272^ 



4 logaV - 24.36875 (1 - a)( - logaV + 5.792010 



(81) 



2 2 

-0.12554 • Cf • i74P4 + (7.16687 - 9.17479 (1 - a)) ■ • i — . 

167r^ 167r^ a 
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(a) vertex 



(b) sail 





(c) sail 



(d) operator tadpole 





(e) leg self-energy 
(excluding tadpoles) 



(f) leg self-energy 

(excluding tadpoles) 





(g) leg tadpole 



(h) leg tadpole 



Figure 1: The diagrams needed for the one-loop renormaUzation of the lattice operators. 



Note that gauge invariance forces the terms proportional to (1 — a) to be the same as, for 
example, in the case of Wilson or overlap fcrmions. This is an important check of the cor- 
rectness of our calculations. The two terms proportional to 74/a arising from the tadpole and 
the sunset diagrams do not cancel - they actually reinforce each other. However, the parts 
proportional to (1 — a) cancel exactly, as required by gauge invariance. 



A. 2 Bilinears 

Here wc list the results for the individual vertex diagrams for the scalar density, as well as 
the vector and tensor currents. As a consequence of chiral symmetry, the vertex correction 
for the pseudoscalar density is identical to that of the scalar density. The same is true for the 
vector and axial- vector currents. 

For the scalar and pseudoscalar densities the result for the vertex correction is 



As = 



167r2 



- 4 log a V + 24.36875 (1 - a) ( log a^p^ - 5.79201o) 



(82) 



One can see that the breaking of hyper-cubic symmetry does not induce any mixing. 
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For the vector current the vertex diagram yields 



90 



■ log aV + 10.44610 - 2.88914 • 5^ + (1 - a) ( log aV - 4.792010 



(83) 



This result shows cxphcitly that the spatial and temporal components of the vector (and also 
those of the axial-vector) current receive different radiative corrections. This is a consequence 
of the breaking of hyper-cubic symmetry, and of the special role taken by the temporal di- 
rection. On the other hand, it is encouraging that mixing between the spatial and temporal 
components appears to be absent. Each of these components still renormalizes multiplica- 
tively, and the mixing matrix is diagonal. 

Finally, for the tensor current we obtain the result for the vertex diagram as 



9l 



4.17551 + (1 - a)(logaV - 3.79201o) 



(84) 



Similarly to the scalar and pseudoscalar case, the breaking of hyper-cubic invariance does not 
generate here any extra mixing. It is remarkable that the tensor operator does not appear to 
show any preference for the temporal direction even after (one-loop) renormalization, that is, 
the renormalization constant is the same for each of the six independent components of the 
tensor operator. 
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